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A selfconsistent theory of the current-induced switching of magnetization using nonequilibrium 
Keldysh formalism is developed for a junction of two ferromagnets separated by a nonmagnetic 
spacer in the ballistic limit. It is shown that the spin-transfer torques responsible for current- 
induced switching of magnetization can be calculated from first principles in a steady state when 
the magnetization of the switching magnet is stationary. A steady state is achieved when the spin- 
transfer torque, proportional to bias voltage in the linear response regime, is balanced by the torque 
due to anisotropy fields. The spin-transfer torque is expressed in terms of one-electron surface Green 
functions for the junction cut into two independent parts by a cleavage plane immediately to the left 
and right of the switching magnet. The surface Green functions are calculated using a tight-binding 
Hamiltonian with parameters determined from a fit to an ab initio band structure. This treatment 
yields the spin transfer torques taking into account rigorously contributions from all the parts of 
the junction. The spin-transfer torque has two components, one with the torque vector Tj| in the 
plane containing the magnetizations of the two magnetic layers and another with the torque vector 
Tj_ perpendicular to this plane. It is shown that, in general, Ty and T± may be comparable in 
magnitude and they both tend to finite values independent of the spacer thickness in the limit of a 
thick spacer. T± is shown to be small when the exchange splitting of the majority- and minority-spin 
bands in both ferromagnets tends to infinity or in the case when the junction has a plane of reflection 
symmetry at the center of the spacer. The torques T± and T« are comparable for a Co/Cu/Co(lll) 
junction when the switching Co layer is one or two atomic planes thick. Tj_ is « 27% of T|| even 
for a switching Co magnet of ten atomic planes. Depending on material parameters of the junction, 
the relative sign of T± and Ty can be negative as well as positive. In particular, Tx/Tu < for 
Co/Cu/Co(lll) with switching Co magnet of one atomic plane and Tx/T\\ > for two atomic 
planes of Co. A negative sign of the ratio Tx/Tu has a profound effect on the nature of switching, 
particularly in the realistic case of easy-plane (shape) anisotropy much larger than in-plane uniaxial 
anisotropy. To calculate the hysteresis loops of resistance versus current, and hence to determine 
the critical current for switching, the microscopically calculated spin-transfer torques are used as an 
input into the phenomenological Landau-Lifshitz equation with Gilbert damping. In the absence of 
an applied magnetic field, an ordinary hysteresis loop is the only possible switching scenario when 
Tx/Tu > 0. However, for Tx/Tu < 0, a normal hysteretic switching occurs only at relatively low 
current densities. When the current exceeds a critical value, there are no stable steady states and 
the system thus remains permanently in a time dependent state. This is analogous to the observed 
precession of the switching magnet magnetization caused by a DC current in the presence of an 
applied magnetic field. The present calculations for Co/Cu/Co(lll) show that the critical current 
for switching in the hysteretic regime is ss 10 7 A/ cm 2 , which is in good agreement with experiment. 



INTRODUCTION 

Slonczewski proposed a new method of switching the magnetization direction of a thin film by means of 
a spin-polarized current. The current is spin-polarized by passing through a thick layer of a ferromagnetic 
metal, whose magnetization is assumed to be pinned, subsequently passing through a nonmagnetic metallic 
spacer layer and then through a thin magnetic switching layer into a nonmagnetic lead. Early related 
theoretical work is due to Berger j2| ■ Switching of the magnetization is accompanied by a change in resistance 
(CPP GMR) and the effect has been observed experimentally by studying hysteresis loops in resistance versus 
current plots for pillar systems Jumps in the hysteresis curve occur between steady states of constant 
current and static magnetization, just as in the Stoner-Wohlfarth [J] theory of field-switching jumps occur 
between equilibrium states. We have formulated a first-principle theory of current-induced switching based 
on this idea. As in the Stoner-Wohlfarth theory, we assume that the switching magnet remains single domain 
during the switching process. 
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One of the main aims of this paper is to calculate hysteresis loops of resistance versus current from first 
principles for a much more general situation than has been considered previously. In previous treatments 
a uniaxial anisotropy field was introduced in the switching magnet with its direction parallel to the mag- 
netization of the polarizing magnet 0,0]. I n this case, there are only two steady states in which the 
magnetizations of the polarizing and switching magnets are either parallel or antiparallel. It is for this rea- 
son that the steady-state approach has not previously been further developed. However, in real experiments 
on pillar structures shape anisotropy due to variable shape of pillar cross sections means that the direction of 
the anisotropy field in the switching magnet is not simply related to the direction of the magnetization of the 
polarizing magnet. Furthermore, in some experiments ||J, an external field is also applied so that more 
general orientations of the magnetizations occur in the steady states. It is, therefore, essential to consider 
a completely general case when the uniaxial anisotropy field makes an arbitrary angle 9 with the polarizing 
magnetization. We also include the easy plane anisotropy, which is always large in layered magnets, and 
investigate fully its consequences. 

In a steady state there is a balance between the spin-current torque, acting on the switching magnet due to 
the spin-polarized current, and the torque due to anisotropy and external fields. In our general first-principle 
treatment, two components of the spin-current torque appear naturally, one with the torque vector Zji in 
the plane containing the magnetizations of the two magnetic layers ('in-plane' torque) and another with the 
torque vector Tj_ perpendicular to this plane ('out-of-plane' torque). Slonczewski [lj considered only Tji and 
it is generally believed that Tj_ is negligible. We start the presentation of our results in Sec. 6 by deriving 
from the general Keldysh formalism the results of Slonczewski's original calculations [lj, in which only Ty 
appears. It will be seen that this result is not always valid but is just an artefact of Slonczewski's simple 
model. In fact, we shall show that in some cases T± is dominant and that, even when small, T± is essential 
since its importance is strongly enhanced in the presence of easy-plane anisotropy. 

To calculate hysteresis loops for this general scenario, we need to solve the following problems: (i) calculate 
microscopically both the in-plane and out-of-plane components of the spin-current torque; (ii) determine the 
steady states which form the continuous parts of the hysteresis curve; (iii) investigate the stability of such 
states in order to determine critical currents at which the jumps, and hence switching, occur; (iv) calculate 
the resistance of the layered structure along the steady-state paths. 

Within our unified theory, all this can be done for a general layered system with a fully realistic band 
structure. 

A jump in the hysteresis curve occurs at a critical current when one steady state becomes unstable and 
the system seeks out another stable steady state. This is in analogy with the Stoner-Wohlfarth theory of 
field-switching where one deals with equilibrium states instead of the present noncquilibrium steady states. 
As in that theory we do not concern ourselves with the detailed dynamics of the switching. However, we 
identify in this paper certain cases in which one steady state becomes unstable above a critical current 
but there are no other stable steady states available. Under these circumstances the magnetization of the 
switching layer remains perpetually in a time-dependent state. 

In order to study nonequilibrium steady states we use the Keldysh formalism 0, ^j|> 0] described in 
Sec. 3. As pointed out above, the steady state arises from a balance between spin-current torque and 
anisotropy field torque. Hence it is essential to include anisotropy and/or external fields in the Hamiltonian 
of the system from the outset. It is also necessary to treat correctly the on-site electron-electron interaction 
which is responsible for the spontaneous magnetization of the polarizing and switching magnets. This is 
achieved by insisting that the local exchange field is in the direction of the local magnetization, which 
is the essential feature of self-consistent field approximations such as unrestricted Hartree-Fock (HF) and 
local spin density (LSD A). Such a treatment respects the spin-rotational symmetry of the ferromagnet in 
the absence of external fields. Beyond this we do not need to introduce a self-consistent treatment of the 
Coulomb interaction explicitly, although bulk LSDA calculations underlie the band parameters and exchange 
splittings used in our calculations. In our Keldysh approach the direction of magnetization in each atomic 
plane of the switching magnet is determined self-consistently in the steady state by the requirement that 
the magnetization of a given atomic plane is parallel to the exchange field in that plane. The relationship 
between this approach and the more intuitive one of balancing torques is discussed in Sees. 2 and 5. 

The treatment described above enables us to determine all possible steady states of the system and the 
next step is to investigate their stability. We do this by introducing the spin-current torques, calculated 
microscopically as functions of magnetization direction, and anisotropy torques into a Landau-Lifshitz equa- 
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tion of motion for the magnetization including Gilbert damping. We linearize the equation of motion about 
the steady-state solution to obtain the conditions for stability. 

Finally, we construct hysteresis curves from continuous steady-state paths and jumps at points of 
instability. 



THEORETICAL MODEL 

The layer structure we consider is shown in Fig.l. It consists of a 




Y 



current 

FIG. 1: Schematic picture of a magnetic layer structure for current-induced switching (magnetic layers are darker, non-magnetic 
layers lighter). 

semi-infinite polarizing ferromagnet with magnetization P, a nonmagnetic metallic spacer with N atomic 
planes, a switching magnet with M atomic planes, and a semi-infinite nonmagnetic lead of the same material 
as the spacer. Each layer is described by a tight-binding model, in general multiorbital with s, p, and d or- 
bitals whose one-electron parameters are fitted to first-principle bulk band structure, as discussed previously 
[l2j| . The Hamiltonian is, therefore, of the form 

H = H + Hi n t + H anisi (1) 
where the one-electron hopping term Ho is given by 

H o = Yl 53 T ^^y{k\\)c^ mlia c k ^ nva , (2) 
where ct creates an electron in a Bloch state, with in-plane wave vector fen and spin a, formed from 

K^mna x II r 

a given atomic orbital /i in plane m. rF nt is an on-site interaction between electrons in d orbitals which 
leads to an exchange splitting of the bands in the ferromagnets and is neglected in the spacer and lead. 
The magnetization of the polarizing magnet is assumed to be pinned in the (z,x)-plane, making an angle 9 
with the z axis, as shown in Fig.l. H an i S contains effective fields in the switching magnet corresponding to 
uniaxial H u and easy-plane H p anisotropics. It is given by 

Hants = — ^ ' S n .H a, (3) 

n 

where S n is the operator for the total spin angular momentum of plane n and 



Ha = H u + H 



(4) 
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H u and H p are given by 

H u = (e z . < >)H u0 e z (5) 

H p = -(e y . < S„ >)H pQ e y , (6) 

where < > is a unit vector in the direction of the thermal average of S n , and e x , e y , e z are unit vectors 
in the direction of the axes shown in Fig.l. H u o, H p q measure the strengths of the uniaxial and easy-plane 
anisotropics and have dimensions of frequency. These quantities may be converted to a field in tesla by 
multiplying them by ft/2/ig = 5.69 x 1CP 12 . We assume that anisotropy fields are uniform throughout the 
switching magnet but it would be easy to generalize to include, for example, a surface anisotropy. The spin 
angular momentum operator S n is given by 

Sn = 2^^^ C ^ir^T' C fc||«A'l- )cr ^ Cfc ll™^T' Cfc ll n ^^ T ( 7 ) 

fe||Ai 

and the corresponding operator for spin angular momentum current between planes n — 1 and n is 

3n-i =~2 X! r ( fe ll)"^"- 1 A'( c I| r i/T' C fc||»i'l' ) ' T ( Cfc ll"- 1 A l T' c fc||"-iA'l) T + h - c - ( 8 ) 

fc|| ilv 

Here, cr = {a x , a y , er z ), where the components are Pauli matrices and Eq.(8) yields the charge current 
operator if \<j is replaced by a unit matrix multiplied by the electronic charge e/H, where e is the electronic 
charge (negative). 

All currents flow in the y direction, perpendicular to the layers, and the components of the vector j 
correspond to transport of x, y, and z components of spin. The rate of change of S n in the switching magnet 
is given by 

ihS n = [S n ,H ] + [S n , H ams ] (9) 

since the spin operator commutes with the interaction Hamiltonian Hint- 
It is straightforward to show that 

[S n ,H Q ] = ih(j n _ 1 -j n ) (10) 

and 

[S n ,H anis ]=-ih(H A x S n ), (11) 

In a steady state, the magnetization is time- independent so that < S n >= 0. Hence 

< J«-i > - < in >=H A x < S n > . (12) 

The left-hand side of Eq.(12) corresponds to the rate of transfer of spin angular momentum to plane n in the 
steady state. Thus Eq.(12) shows explicitly how, in the steady state, this spin-transfer torque is balanced by 
the torque due to anisotropy fields. The concept of spin-transfer torque was first introduced by Slonczewski 

III 

KELDYSH FORMALISM 

In this section we show how to calculate the spin current < j n _x > and spin density < S n > in the non- 
equilibrium steady state and verify that they are related by Eq.(12). To produce a spin-polarized current in 
the system we apply a bias Vb between the polarizing magnet and the lead. To use the Keldysh formalism 
0, [ljj, to calculate < J n _i > and < S n > we consider an initial state at time t — — oo in which the 
hopping integral T n ^, n -ifi between planes n — 1 and n is switched off. Then both sides of the system are in 
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equilibrium but with different chemical potentials \il on the left and hr on the right, where Mfl = e ^fe- 
The interplane hopping is then turned on adiabatically and the system evolves to a steady state. The cleavage 
plane, across which the hopping is initially switched off, may be taken in either the spacer or switching layer 
or in the lead. Fig.l shows the situation when the cleavage plane is between atomic planes n — 1 and n in 
the switching magnet. In principle, the Keldysh method is valid for arbitrary bias Vb but here we restrict 
ourselves to small bias corresponding to linear response. This is reasonable since for larger bias electrons 
would be injected into the switching magnet far above the Fermi level and many-body processes neglected 
here would be important. Furthermore, in metallic systems the bias will never be large. 
Following Keldysh , 01 > we define a two-time matrix 

G+ L (t,t')=i<cl(t')c R (t)>, (13) 

where R = (n, v, a 1 ) and L = (n — cr), and we suppress the fcy label. The thermal average in Eq.(13) is 
calculated for the steady state of the coupled system. The matrix G~^ L has dimensions 2m x 2m where m is 
the number of orbitals on each atomic site, and is written so that the m x m upper diagonal block contains 
matrix elements between | spin orbitals and the m x m lower diagonal block relates to J, spin. 2m x 2m 
hopping matrices Tlr and Tut are written similarly and in this case only the diagonal blocks are nonzero. 
If we denote Tlr by T, then Tel — . We also generalize the definition of cr so that its components are 
now direct products of the 2x2 Pauli matrices a x ,a y , a z and the m x m unit matrix. The thermal average 
of the spin current operator, given by Eq.(8), may now be expressed as 

< J'n-i >=lj2 Tr ^ G RL(t,t)T-G+ R (t,t)T^] ( r}. (14) 
fey 

Introducing the Fourier transform G + (w) of G + (t, t'), which is a function of t — t', we have 

< 3n-x >= \ [ ^ T r{[G + MT (15) 
fe|. 

Again, the charge current is given by Eq.(15) with icr replaced by the unit matrix multiplied by e/H. 

Similarly, the total spin angular momentum on atomic planes on either side of the cleavage plane, in the 
non-equilibrium state, is given by 

< Sn-i >= -\^T. J ^Tr{G + LL (cj)a} (16) 



< S n >= ~itiJ2 J ^Tr{G+ R {u)*}. (17) 
fe|. 

Following Keldysh [ljj, we now write 

G + AB ^) = \{F AB + G a AB -G\ B ), (18) 
where the suffices A and B are either R or L. Fab(u) is the Fourier transform of 

F AB (t,t') =-i< [c A {t),c B {t')]- > (19) 

and G°, G r are the usual advanced and retarded Green functions Note that in [j|and 0] the definitions 
of G° and G r are interchanged and that in the Green function matrix defined by these authors G + and G~ 
should be interchanged. 

Charge and spin current, and spin density, are related by Eqs.(15)-(17) to the quantities G a , G r , and F AB . 
The latter are calculated for the coupled system by starting with decoupled left and right systems, each in 
equilibrium, and turning on the hopping between planes L and R as a perturbation. Hence, we express G a , 



G 



G r , and Fab in terms of retarded surface Green functions <?l — 9ll, 9r = 9rr for the decoupled equilibrium 
system. The final result for the spin angular momentum on plane n to the right of the cleavage plane is 

< S n >=< S n >i + < S n > 2 , (20) 

where the two contributions to the spin angular momentum < S n >i and < S n >2 are given by 

< 5„>i= J duJ ImTr{Ag R a}[f(to - n L ) + f(u - » R )} (21) 

fen 



< S n > 2 = - 




Here, A = [1 - gj^T^ g^T]^ 1 , B = [1 - g^Tt^T]" 1 , and /(w - m) is the Fermi function with chemical 
potential ^ and (1l — M-R = e ^b- To obtain < S„_i >, defined by Eq.(16), we must interchange L and i?, and 
T and 7^, everywhere in Eq.(20)-(22). In the linear-response case of small bias which we are considering, 
the Fermi functions in Eq.(22) are expanded to first order in V h . Hence the energy integral is avoided, being 
equivalent to multiplying the integrand by eVj, and evaluating it at the common zero-bias chemical potential 
Mo- 

As shown in Fig.l, the magnetization P of the polarizing ferromagnet is assumed to be fixed in the (z, x) 
plane and makes an angle 9 with the z axis, which is the direction of the uniaxial anisotropy field in the 
switching magnet. When a bias is applied, spin-polarized current flows through the switching magnet and 
exerts a torque on its magnetization. This torque is in competition with the torque due to the anisotropy 
field and causes the spin < S n > in a given atomic plane n to deviate from the anisotropy axis. In the steady 
state < S n > settles in a definite direction specified by the angles a n , <fi n shown in Fig.l. To determine these 
angles, we assume the exchange field A„ in plane n is in the direction (a n , <fr n ) and apply the self-consistency 
condition 

A„x < >=0. (23) 

This condition guarantees that the local magnetization is in the direction of the exchange field, as it should be 
in the unrestricted Hartree-Fock approximation mentioned in Sec.l. As with anisotropy fields, the exchange 
field A n is defined as an angular frequency so that fiA n is the energy to reverse the spin on plane n. More 
precisely, the spin-dependent part of the on-site energy on plane n is given by — {l/2)h(H a + A„).<r. We 
assume that | A„ | always takes its bulk value. 

Following the method outlined for obtaining Eq.(20), similar expressions in terms of retarded surface Green 
functions may be obtained for the spin currents < j n _ 1 > and < j n >. Writing again < j n >=< j n >i 
+ < j n >2, we obtain 




< J„_i >2= ^ Y,j duReTr{[g L TABg R Tt - AB + ±(A + B)]tr}[f(u - /i L ) - f(u> - ^ R )\. (25) 
fey 

By considering the changes in g^, g R when the cleavage plane is moved one atomic plane to the right, it is 
straightforward to show that 

< 3n-\ >~ <3n >= ( H A + A„)x < S n > . (26) 

This equation holds for a steady state with arbitrary exchange fields A„ which do not necessarily satisfy 
the self-consistency condition (23). When the self-consistency condition (23) is satisfied, we recover the 
steady-state result 



< Jn-l > - < Jn >= ^ X < S * > 



(27) 
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which was derived earlier [Eq.(12)] purely from considerations of the spin-rotational symmetry of the electron- 
electron interactions. This verifies the consistency of the Keldysh formalism combined with the unrestricted 
Hartree-Fock approximation. 

It follows from Eqs.(26) and (27) that all components of spin current are conserved within the spacer and 
lead, where Ha = 0, A n = 0, with or without self-consistency. Furthermore, it follows from Eq.(27) that in 
the self-consistent steady state the component of spin current in the direction of the anisotropy field Ha is 
conserved throughout the system, as of course is the charge current. If Eq.(27) is summed over all planes in 
the switching magnet, we obtain 

< 3 spacer >~< 3 lead >= H A X < S t ot >, (28) 

where < j spacer >, < 3i ea d > are the spin currents in the spacer and lead, respectively, and < Stot > is the 
total spin angular momentum of the switching magnet. This shows how the total spin transfer torque acting 
on the switching magnet is balanced by the torque exerted by the anisotropy field on the total moment. 

We have separated in Eq.(20) the spin angular momentum < S n > into two parts < S n >i and < S n >2- 
It is clear that < S n >2 is proportional to the applied bias Vb to the first order and for zero bias — (J-r) 
only < S n >i remains. The spin transfer torque < j n > — < j n _ ± > similarly splits into two parts 
(Eqs.(24),(25)) in such a way that Eq.(26) holds for each component separately: 

< j n -i >i~ < 3„ >i= [Ha + A„)x < S n >i, i = 1,2. (29) 

Only the first part < j n _j >i — < j n >i is nonzero at zero bias. It corresponds to spin currents which 
mediate exchange coupling, either between the two magnets across the spacer or between atomic planes in 
the switching magnet. Consequently, at zero bias the spin current in the lead is zero. It is easy to verify that 
the expressions for interlayer exchange coupling^ derived here, using the Keldysh formalism, agree precisely 
with those obtained earlier by other methods 

The results of this section show the great advantage of the Keldysh formalism, even within the linear 
response regime. Spin currents at zero bias, corresponding to exchange coupling, transport spin and particle 
currents, and spin densities are all calculated in a unified way. Relationships between these quantities, such 
as Eq. (29), are then easily derived. In the standard linear response theory of Kubo zero-bias quantities 
cannot be calculated and different response functions would have to be introduced for calculating currents 
and spin density response at finite bias. 



APPLICATION TO A SWITCHING MONOLAYER 



In the general theory outlined in Sec. 3 the steady-state spin orientation of each atomic plane n of the 
switching magnet must be determined self-consistently. In this section we first consider the simplest case 
of a single orbital on each site and when the switching magnet is a single atomic plane. In this case there 
is no interplane exchange coupling in the switching magnet to consider and we assume that the spacer is 
sufficiently thick for the zero-bias exchange coupling between the two ferromagnets to be negligible. For a 
given bias V&, the direction (ao, </>o) of the steady-state orientation of the switching magnet moment < S > 
is determined self-consistently from Eq.(23) with the cleavage plane immediately to the left of the switching 
plane so that < S n >=< S >. It is convenient to determine (a ,<j) ) in two steps. The first step locates a 
'universal path' on a unit sphere, independent of Vb, on which the self-consistent solutions for any given V 
must lie. In the second step the bias Vb required to stabilize the magnetization in a given direction (ao,</>o) 
is determined as a function of ao, say. To establish this program, we write Eq.(23) as 

(H A + A)x < S >= H A x < S >, (30) 

where A is the exchange field of the switching layer in the direction (ao, 4>q). Splitting < S > into two parts 
as in Eq.(20), this becomes 

{H A + A)x<S>! +(H a + A)x < S > 2 = H A x < S > . (31) 

Hence using Eq.(29) we have 

< Jspacer >1 ~ < 3 lead >1 + < jspacer >2 - < 3 lead >2= H A x < S > . (32) 
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The first two terms on the left correspond to exchange coupling torque which, as discussed above, is assumed 
to be negligible compared with the anisotropy torque. This is justified for thick spacers since the interlayer 
exchange coupling tends to zero as the spacer thickness tends to infinity. The last two terms on the left 
correspond to the spin transfer torque T, which is proportional to bias Vb, and the right-hand side of Eq.(32) 
is — T a-, where T a is the torque exerted by the anisotropy field on the switching magnet. We shall see that, 
in contrast to the exchange coupling torque, T remains finite as spacer thickness tends to infinity. Hence 



T=< j 



spacer 



>2 - < 3 lead >2= H A X <S> = —T/ 



and, in particular, 



T.H, 



0. 



(33) 



(34) 



The bias Vb now cancels and this equation determines the universal path described above. Eq.(33) determines 
the bias required to stabilize any particular point on this path of possible steady states. 

We conclude this section with one example in which, for simplicity, we retain only uniaxial anisotropy, 
this field being chosen in the z direction. We use a single-orbital tight-binding model whose lattice is taken 
to be simple cubic with layering in the (010) direction. The nearest-neighbor hopping parameter t is taken 
to be the same throughout the system. The on-site energy in the spacer and lead is taken as V sp , the zero of 
energy being at the common Fermi level for zero bias. In this example the on-site energy V sp is also taken for 
majority spin in the ferromagnets (perfect matching in the majority-spin channel). The on-site energy for 
minority spin in the ferromagnets is taken as V sp + fiA, where fiA is the exchange splitting. The matching 
of spacer and majority spin bands is similar to the situation in Co/Cu. We take V sp = 2.3, fiA — 0.7 in 
units of 2t. Furthermore, we take the uniaxial field parameter H u q — 1.86 x 10 10 sec _1 which corresponds 
to a field of 0.106T. We also take a general value 9 — 2 radians of the angle between the polarizing magnet 
moment and the direction of the uniaxial anisotropy axis of the switching magnet. To determine the torque 
T which appears in Eq.(34) for the universal path, we need to calculate the Green functions ql, gn which are 
required in Eq.(25). This is done by standard adlayering methods described previously [12j. At this stage, 
the anisotropy field is included in the calculation of all the Green functions. Fig. 2 shows the calculated 
universal path of 4> versus a for the specified parameters. The bias Vb required to yield a steady-state 




FIG. 2: Universal path of <j> versus a. 



magnetization at a given point of this path is plotted as a function of a in Fig. 3, where we have assumed the 
band width 12t — 6eV. Positive bias corresponds to a drop in voltage between the polarizing magnet and 
the lead. The correspondence between the two curves in Fig. 2 and Fig. 3 is indicated by letters P and Q. 
The discussion of stability of these steady states and the interpretation of Fig. 3 is postponed to section 7. 
The method of calculating steady states used in this section becomes more complicated when the switching 
magnet contains several atomic planes since the moments < S n > of all planes must be determined self- 
consistently. This entails the inclusion of the exchange stiffness between atomic planes of the switching 
magnet which is contained in Eq.(29) with i = 1. To address this problem, we introduce in the next section 
the simplifications required to derive from first principles the convenient 'standard model' used by previous 
authors Q, 0, H- 



angle a 



FIG. 3: Bias Vf, required to stabilize the switching magnet moment at an angle a on the universal path. 

THE STANDARD MODEL 

In the Keldysh method used above it is essential to include the anisotropy field Ha in the Hamiltonian 
to obtain a non-trivial steady state. In the absence of Ha it follows from Eq.(27) that in the steady state 
all components of spin current are conserved everywhere so that there are no spin-transfer torques. Hence 
the only steady state is the trivial one in which the switching magnet is aligned parallel or antiparallel to 
the polarizing magnet. Previous authors 0, § did not consider a steady state but calculated spin-transfer 
torque as a one-electron problem with the exchange fields of the polarizing and switching magnets at a fixed 
assumed angle. In a second independent step, these authors [lj, |5j balance the spin-transfer torque against 
the torque due to anisotropy field in the context of a Landau-Lifshitz equation. This approach is what we 
call the standard model (SM). In this section we show how to arrive at the SM by making some simplifying 
approximations in our self-consistent approach. 

We begin with the monolayer switching magnet of Sec. 4. In Eq.(33) the spin-transfer torque T is calculated 
in the presence of Ha and the spin < S > is the self-consistent moment. To obtain the SM we must neglect 
Ha in the calculation of T and replace < S > by its nonself-consistent value in the direction of the assumed 
exchange field of the switching magnet. These approximations are both reasonable provided the exchange 
field is much stronger than the anisotropy field, which is satisfied for a ferromagnet such as Co. This follows 
since the Green functions which determine both T and < S > depend on the total field A + Ha- Thus in 
the SM the spin-transfer torque is calculated as a function of the angle between the magnetizations without 
solving the self-consistency problem. Furthermore, equating it to the anisotropy torque as in Eq.(33) is 
equivalent to calculating a steady state of the Landau-Lifshitz equation. The justification of the SM for a 
switching magnet with more than one atomic plane is more subtle. 

The self-consistency condition (23) must be satisfied for each plane in the switching magnet. It may be 
written [c.f. Eq.(30) for the monolayer] 

(H A + A n )x <S n >=H A x <S n > (35) 

and, using again < S n >=< S n >i + < S n >i and Eq.(29), we obtain 

< 3n-l >1 - < 3n >1 + < 3n-l >2 ~ < j n >2 = H A X<S n >. (36) 

The first two terms contain the interlayer exchange coupling, which is neglected as in the monolayer case, 
and interplane exchange coupling within the switching magnet. To clarify the argument, we write notionally 
this last contribution in terms of local exchange stiffness D n between atomic planes n— 1 and n. Hence from 
Eq.(36) 

< 3n-l >2 - < j n >2= H A x < S n > +D n _ 1 (S n - 1 x S n ) - D n (S n x S n+1 ). (37) 
On summing over all planes n in the switching magnet the internal exchange coupling torques cancel and 
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we have 

< 3 spacer >2~ < 3 lead >2= H A X ^ < S n > . (38) 

n 

In the fully self-consistent solution of Eq.(35) the exchange field A„ is parallel to the local moment < S n > 
but the A n are not collinear. To proceed to the SM we must assume that all the A„ used to calculate the 
spin-transfer torque on the left of Eq.(38) are equal, to A, say. Furthermore, we assume, as in the case of 
the monolayer, that we can neglect Ha in the calculation of the spin-transfer torque and that all < S n > 
are in the direction A with magnitude equal to the ground state moment. In making the approximation 
A„ w A we have failed to satisfy Eqs.(37) individually as is required for full self-consistency. 

To show that this is not serious for a ferromagnet such as Co, consider the following argument. If we 
use the uniform value of A determined from Eq.(38), as described, to calculate the spin-transfer torque in 
Eq.(37) and assume all < S n > are in the direction of A, the last two terms of Eq.(37) are zero and the 
equations are far from satisfied. However, since the exchange stiffness constants D n of a ferromagnet such 
as Co are large one need only introduce small deviations of A„ from the uniform A, and consequently small 
deviations of < S n > from uniformity, to satisfy the self-consistent equations (37). This is true because 
the spin-transfer and anisotropy torques are insensitive to these small deviations. The ability of the SM to 
simulate the fully self-consistent solution accurately has been verified numerically for a switching magnet 
with two atomic planes using the single-orbital model of Sec. 4. 



TWO COMPONENTS OF THE SPIN-TRANSFER TORQUE IN THE STANDARD MODEL 

In the calculation of the spin-transfer torque T within the standard model the anisotropy field is neglected 
so that T depends only on the angle tp between the magnetization P of the polarizing magnet and the 
assumed exchange field A of the switching magnet. As in Fig.l, the magnetization P of the polarizing 
magnet is in the (z,x) plane, making an angle 6 with the z axis, and, for convenience, we choose the exchange 
field of the switching magnet to be in the z direction so that ip = 9. The torque T in the SM is given by 

T =< 3 spacer >2 ~ < 3 lead >2 (39) 

where the right-hand side is related to the total bias-induced spin S2 = < S'n >2 by Eq.(29), summed 
over n, with Ha neglected and A„ = A. It follows that 

T = A x < S n >2= A X S 2 . (40) 

n 

Here, < S n >2 is given by Eq.(22) with Ha neglected in the Green functions <?£, g^. The three components 
of S2 = {S2X, S2y, S2z) are related to the three Pauli matrices a x ,a y ,a z in Eq.(22). Clearly, from Eq.(40) 
the z component of torque is zero so that we can write 

T=(T x ,T y ,0). (41) 

The 'in-plane' component Tj| = T X1 where 'in-plane' refers to the (z,x) plane containing P and A, is given 
by Tj| = — AS2?/ an d the 'out-of-plane' component T± = T y is given by T± = AS^. The quantities 
S2X and S2 y represent small deviations of the switching magnet moment from the direction of its exchange 
field. These spin components are referred to by previous authors 0|> 01 as ' s P m accumulation'. In the 
self-consistent steady-state treatment of Sec. 4 such deviations do not occur because the exchange field is 
always in the direction of the local moment. In our view, time-independent spin accumulation Six, S2y in 
the ferromagnet is a non-physical concept which, however, we may define formally as the ratio of torques Tj| , 
T±, to the exchange field A. It is remarkable that, as shown in Sec. 5, the SM in which this concept arises 
provides a convenient and frequently accurate method for calculating spin-transfer torque. Time-dependent 
spin accumulation in the ferromagnet in a non-steady state could be a valid concept. However, these time- 
dependent spins would produce time-dependent exchange fields which would excite the whole spin system. 
This would require a many-body treatment going beyond the unrestricted Hartree-Fock approximation which 
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is adequate for the steady state. Spin accumulation < S x >, < S v >, < S z > proportional to the bias exists 
in the spacer even in the steady state and has real physical significance. Calculation of this effect will be 
published in a succeeding paper. 

The spin transfer torque can be calculated either directly from Eq.(39) or from Eq.(40). However, the 
latter would require calculating < S n >2 for each atomic plane of the switching magnet so that the direct 
method is obviously preferable. 

We begin with an exactly solvable one-band model which we can connect with previous work yj. This 
model is related to the one described at the end of Sec. 4, where the switching magnet is a single atomic plane 
and there is perfect matching between the spacer band and the majority spin bands in both ferromagnets. 
To obtain analytical results in this first example, we also assume that the exchange splitting A — > oo both in 
the polarizing and switching magnets. In fact, once the bottom of the minority-spin band is well above the 
Fermi level, the results are rather insensitive to the magnitude of A. Such a system is sometimes referred 
to as a half- metallic ferromagnet and is the first case considered by Slonczewski in his original paper In 
the limit A — ► oo the SM model is exact since the moment of the switching magnet cannot deviate from 
the exchange field and the self-consistency condition (23) is automatically satisfied. We, therefore, calculate 
the spin-transfer torque in the absence of anisotropy field. Clearly, owing to the infinite exchange splitting, 
the only spin current in the lead corresponds to the z component of spin and the z-spin current is equal to 
the charge current (multiplied by H/2e). It turns out in this model that the y-spin current in the spacer, 
which is equal to the torque Tj_ since the corresponding current in the lead is zero, vanishes. Thus only 
Slonczewski torque Tm survives and is given by 

T = -p — V f2 ^°~ go) 2 sin t (AO) 

11 x 8n ^ \1~ t 2 g (a + bcosip) | 2 ' 1 1 

fc,l 

Here, go = go(fc||,0) is the majority spin surface Green function for the semi-infinite ferromagnet, or equiv- 
alently for the semi-infinite ferromagnet with an overlayer of the matching spacer. The Green function go is 
evaluated at energy u — 0, the common Fermi level of the unbiased system. Also, t is the hopping parameter 
introduced in Sec. 4 and a and b are given by 

a=^{g + A N ) ; b=^(g -A N ), (43) 

where An = sin Nk±d/[t sin(iV + l)k±d] with 

k ± d = cos" 1 [(y s T +2t(cosk x d + cosk y d))/2t] (44) 

and feii = (k x , k y , 0); d is the interatomic distance, N is the number of atomic planes in the spacer, and Vj 
is the on-site potential in the majority-spin band of the switching magnet. The corresponding expression for 
the charge current is very similar and we find that 

Tii = — — - tan(-0/2) x (charge current). (45) 
2 | e | 

This is precisely the Slonczewski result for the analogous parabolic band model. It should be noted that 
the torque Tit goes to zero for ip — > n since the charge current for a halfmetallic magnet contains a factor 
1 + cost/i. 

The interesting result that T± = for this model may be traced to an effective reflection symmetry of 
the system about a plane at the center of the spacer. Although the present system appears asymmetric the 
infinite exchange splitting makes it equivalent to a symmetric system with semi-infinite switching magnet. 
More generally, we find, certainly for a one-band model with arbitrary parameters, that the y-spin component 
of the spin current in the spacer always vanishes for a system with reflection symmetry. In general, however, 
the y-spin current in the lead is non-zero so that Tj_ ^= even for a symmetric system. The result Tj_ = 
for the above model is, therefore, a very special one due to the artefact of a very large exchange splitting in 
the ferromagnets. 

In the second set of examples we consider several cases, within the one-band model, where the exchange 
splittings in the ferromagnets are finite. Simple formulas, such as Eq. (42), for the torques are no longer 
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TABLE I: Parameters for one-band models. VJ , V£ are on-site potentials for majority and minority spin in the 
polarizing magnet; Vj , V} are on-site potentials for majority and minority spin in the switching magnet; V sp is the 
on-site potential in the spacer and the lead; N, M are the numbers of atomic planes in the spacer and switching 
magnet, respectively. 
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available and they must be calculated numerically. In all the examples the calculated torques are per surface 
atom. In all cases, we retain the geometry of the first example. Table I lists the parameters for all the cases 
considered. 

All potentials in Table 1 are in units of 2t and the Fermi energy = 0. 

Figure 4 (a) shows the calculated torques Tji per surface atom (in units of eVj,) as a function of the angle 
i\) for the models with parameter sets (a)-(d) of Table 1. In case (a) T\\ dominates, as expected from the 
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FIG. 4: Dependence of the spin-transfer torques Ty (a) and T±_ (b) on the angle ip for the model with parameter sets (a)-(d) 
of Table 1. The torques are in units of eVJ,. 

large exchange splitting Vj — V} in the switching magnet which approaches the infinite exchange splitting 
of our first exactly solvable model. The angular dependence of both torques is clearly dominated by a sin?/) 
factor (c.f. Eq.(42)) although some distortions are apparent. In case (b) Ty and T± are of almost equal 
strength. This is case where the two ferromagnets are of the same material and the bottom of the minority 
spin band is exactly at the Fermi level. This simulates well the situation in Co/Cu and we shall see presently 
that in realistic calculations for this system the torques Ty and T± are again similar in magnitude. The 
parameters of case (b) were used previously in Sec. 4 as an example of a fully self-consistent calculation of 
steady states. In case (c) T± is larger than T». It is interesting that this occurs for smaller exchange 
splitting in the switching ferromagnet. In cases (b) and (c) the angular dependence of the torques is hardly 
distorted from the sin?/> form. In cases (a), (b), (c) the two torques Ty, T± have the same sign. In case (d) 
they have opposite sign and almost equal magnitude. In examples (a)-(d) the switching magnet consists of 
one atomic plane. In case (e) shown in Fig. 5 we use the same parameters as in case (c) but the number of 
atomic planes in the switching magnet varies between 1 and 10 and the angle tp = n/2. It can be seen from 
Fig. 5 that, contrary to popular belief, the out-of-plane torque Xj_ dominates over Ty for small thicknesses 
of the switching magnet and remains 50% of Ty at M = 10 at. planes. For all thicknesses of the switching 
magnet, Tj| and T± have the same sign. We have already seen in case (d) that this is not always the case. 
So far we have kept the number of atomic planes in the spacer at 20 but we must now highlight an 
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FIG. 5: Dependence of the spin-transfer torque Tji and T± on the thickness of the switching magnet M for ip = 7r/2 and the 
parameter set (c) of Table 1. The torques are in units of eVj,. 



important and surprising result concerning the dependence of torque on N . In Fig. 6 we show the torques 
for the parameter set (b) and if> = tt/2 plotted as functions of N. It is clear that they both oscillate but 
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FIG. 6: Dependence of the spin-transfer torque Tn and Tj_ on the thickness of the spacer N for ip = n/2 and the parameter 
set (b) of Table 1. The torques are in units of eVf,. 

tend to constant values as N — ► oo in our ballistic limit. In zero bias all spins of the system lie in the 
(z,x) plane and it is, therefore, not surprising that when charge current flows in nonzero bias there are finite 
z-spin and x-spin currents for arbitrary spacer thickness. The constant value of y-spin current as N — > oo, 
this being associated with T± , may seem more surprising since there is no y-spin density in zero bias. Of 
course, for -0 ^ 0, there is a y-spin current even in zero bias corresponding to interlayer exchange coupling, 
but this effect is not associated with charge transport and decays as 1/N 2 with increasing spacer thickness. 
The relations Tn = — AS^y, Tj_ = AS^ derived in the beginning of this section within the standard model 
show that bias- induced in-plane spin density is related to out-of-plane torque and vice versa. It is, therefore, 
inevitable that both torques will exist. 

As a final example of these spin-transfer torque calculations we consider a fully realistic multi-orbital 
model of fee Co/Cu(lll) with tight-binding parameters fitted to the results of first-principles band structure 
calculations, as described previously Referring to Fig.l, the polarizing magnet is a semi-infinite slab 
of Co, the spacer is 20 atomic planes of Cu, the switching magnet contains M atomic planes of Co with 
M—l and 2, and the lead is semi-infinite Cu. Figure 7(a), (b) shows the angular dependences of Tn, Tj_ for 
the cases M — 1 and M = 2, respectively. For the monolayer switching magnet, the torques T± and Tn 
are equal in magnitude and they have the opposite sign. However, for M = 2, the torques have the same 
sign and T± is somewhat smaller than Tn . A negative sign of the ratio of the two torque components has 
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FIG. 7: Dependence of the spin-transfer torque Ty and Tx for Co/Cu/Co(lll) on the angle -0. The torques per surface atom 
are in units of eVt- Figure (a) is for M=l, and (b) for M=2 monolayers of Co in the switching magnet. 

important and unexpected consequences for hysteresis loops as discussed in the next section. Finally, we 
show in Fig. 8 the dependence of T± and Tn on the thickness of the Co switching magnet. It can be seen that 
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FIG. 8: Dependence of the spin-transfer torque Ty and T± for Co/Cu/Co(lll) on the thickness of the switching magnet M 
for tp = tt/3. The torques are in units of eVt- 

the out-of-plane torque T± becomes smaller than Ty for thicker switching magnets. This is the expected 
behavior since our polarizing magnet is semi-infinite Co, so that as the switching Co magnet becomes thicker 
we approach the limit of a symmetric junction for which the y-componcnt of the spin current vanishes and 
the corresponding component in the lead is usually small. However, T± is by no means negligible (27% 
of Ty) even for a typical experimental thickness of the switching Co layer of ten atomic planes. It is also 
interesting that beyond the monolayer thickness, the ratio of the two torques is positive with the exception 
of M = 4. 



STABILITY OF STEADY STATES AND HYSTERESIS LOOPS 



In Sec. 4 we calculated the steady-state orientation of the magnetization of the switching magnet fully self- 
consistently as a function of bias V h (see Fig. 3). These results allow us to calculate the continuous portions 
of the hysteresis loops of resistance versus bias, but to determine where jumps occur we must investigate the 
stability of the steady states. This cannot be done within the standard Keldysh formalism since the dynamics 
of the system, including damping, lies outside its scope. We, therefore, map the microscopic problem onto a 
phcnomcnological time-dependent Landau-Lifshitz (LL) equation with Gilbert damping. An approach based 
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on LL equation has been used previously [l|, [5J but our treatment differs in several important ways. Firstly, 
we use as an input microscopically determined torques Tj|, T±. In previous work only Tj| was considered 
and its magnitude was an adjustable parameter. Secondly, the importance of the steady state concept in 
the context of the LL equation has not previously been fully recognized. Finally, the consequences of easy 
plane anisotropy have hardly been explored; we find that, in fact, it can lead to completely new switching 
scenarios when Xji and T± have the opposite sign. Even in the absence of easy plane anisotropy, Tji and 
Tj_ with the opposite sign may lead to qualitatively different types of switching. 

In Sec. 6 the unit vector in the direction of the switching magnet moment m was always taken in the z 
direction but to discuss the LL equation we must consider m in a general direction. The total spin-transfer 
torque T may be written quite generally as the sum of the two components in the directions of the vectors 
m x p and m x (p x m) , where p is a unit vector in the direction of the magnetization of the polarizing 
magnet. Thus 

T = T ± +T h (46) 

where 

T± =g±{mx p) (47) 



Ty = gn m x (p x m). (48) 

The modulus of both vector products in Eqs.(47), (48) is equal to smip, where ip (0 < ip < ir) is the 
angle between p and m. Since in the SM T± and Ty depend only on the angle ip, it follows that the 
coefficients gj_, g\\ are functions only of ip. To determine g±(ip), g\\{ip), we return to the geometry of Sec. 6 
with p = (sin-0, 0, cosip) and m = (0, 0, 1). It follows that 

Ty =. 9|| W (sin ^,0,0) (49) 



T ± =g ± (ip)(0,sm*p,a). (50) 

Thus the scalar quantities calculated in Sec. 6 are Tji = g\\ sin-0, T± = gj_ simp. Hence the magnitudes and 
signs of g\\(ip), g±(ip) are determined. It is seen in Fig. 4 that the sin-0 factor accounts for most of the 
angular dependence of Tji , Tj_ so that to a good approximation <7|| , g± are constants, proportional to bias. 
The LL equation takes the form 

dm dm 

— — +7mx— — = T (51) 
dt ' dt y ' 

with the reduced total torque T given by 

T = (-H A x < S tot > +T ± + Tj,)/ |< S tot >|, (52) 

where < Stot > is the total spin angular momentum of the switching magnet and 7 is the Gilbert damping 
parameter. Following Sun 5], Eq.(51) may be written more conveniently as 

(l + 7 2 )^=r- 7 mxr. (53) 

We first consider steady-state solutions of this equation but shall return to the full time-dependent equation 
when discussing stability of these states. In the steady state Eq(51) reduces to 

- H A x < Stot > +T± + Ty = (54) 

which is equivalent to Eq.(33). 

The magnetization unit vector p of the polarizing magnet is given by (sin 9, 0, cos 9) and in the phenomeno- 
logical treatment, based on the SM, the magnetization of the switching magnet is uniform in the direction 
m = (sin a cos (p, sin a sin 0, cos a) (see Fig.l). The procedure for finding steady states is exactly analogous to 
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that described in Sec. 4 for the microscopic approach. Thus the universal path on the unit sphere consisting 
of points (a,4>) which correspond to possible steady states, independent of bias, is again given by Eq.(34). 
For given 9 the torque T = T± + Tji is now defined explicitly as a function of a and cf> by Eqs.(47), (48), the 
bias factor in the constants g±, g\\ cancelling as before. Similarly, using Eq.(54), we can plot a against bias 
Vb for the actual steady state in the given bias as in Fig. 3. 

We now return to Eq.(51) to discuss stability of the steady states. The torque T is given by 

r = H u o{(m.e z )m x e z — h p (m.e y )m x e y + vm x (p x m) + vrm x p}, (55) 

where the relative strength of the easy plane anisotropy h p — H p q/H u q, using the notation of Eqs.(5), (6). 
The last two terms correspond to Tji, with strength parameter v, and T±, with strength parameter rv. 
Clearly the reduced bias v is proportional to the actual bias Vb and inversely proportional to the number of 
atomic planes in the switching magnet. Comparing Eqs.(47), (48), and (55) we see that r = g±/g\\- Thus if 
the scalar torques Tj| , T± defined after Eq.(50) have the same sign it follows that r is positive. Conversely, r is 
negative if the torques have opposite sign. We now linearize Eq.(51) about a steady-state solution m = mo, 
which satisfies T = 0, using the local coordinate axes shown in Fig. 9. Thus 
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FIG. 9: Local coordinate axes for the deviation of the magnetization m from its steady-state orientation mo- 



rn = m + £e a + r\e^ (56) 
and the linearized Eq.(53) may be written in the form 

^ = H + Br,; ^ = C£ + Br), (57) 

where A, B, C\ D are functions of o.q, tpo, 6 and the parameters h p , v, r, and 7. Following Sun we have 
introduced the natural dimensionless time variable r = tH u0 /(l + 7 2 ). The conditions for the steady state 
to be stable are 

F = A + D < 0; G = AD-BC>0 (58) 

excluding F = G = 0[l7j. The damping parameter 7 appears in G only as a factor 1 + j 2 and may be 
cancelled in the condition for stability G > 0. This condition becomes 

Q 2 v 2 + (Qvr + cos2a )(Qvr + cos 2 ao) + h p {Qvr(l — 3sin 2 <f) sin 2 a ) + cos2ao(l — 2sin 2 ao sin 2 </>o)} — 
h 2 s\n 2 a sin 2 <p (l - 2sin 2 O sin 2 a ) > 0, (59) 

where Q = p.rriQ = sin 9 sin ao cos 4>o + cosf?cosao- Similarly, the condition F < becomes 

- 2v(l + jr)Q - 7(cos 2a + cos 2 a ) - ^h p (l ~ 3sin 2 4> sin 2 a ) < 0. (60) 

A number of general conclusions can be drawn from these inequalities. However, we first consider the special 
case when the magnetization of the polarizing magnet is in the direction of the uniaxial anisotropy axis of 
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the switching magnet (6 = 0). In this case the equation r = 0, with T given by Eq. (55), shows immediately 
that possible steady states are given by ao = 0, 7r, corresponding to the switching layer moment along the 
axis of the uniaxial anisotropy. These are the only solutions when h p = 0. However, in the presence of 
easy-plane anisotropy (h p ^ 0) there are additional steady state directions of the switching layer moment 
given by sin20o = 2v/h p , cosao = —vr cos 4>q / (cos 4>q + v sin (f>o) . In the pure Slonczewski case of T± = 
(r = 0) it follows that ao = tt/2. The stability conditions (59) and (60) then reduce to cos2</>o < and 
h p (l — 3 sin 2 4>q) > 1 respectively. For practical biases v/h p <C 1 and the solution for 4>q satisfying the first 
stability condition is <fio « ir/2 — v/h p . The second stability condition is then not satisfied. Thus for 8 = 0, 
r = the only stable steady states of interest are given by ao = or it. It is easily seen that in this case 
the inequality (59) is always satisfied and the system becomes unstable when the left-hand side of Eq.(60) 
changes sign at v c = —(1 + \h p )^i, corresponding to Sun's [5j result in zero external field. This shows clearly 
that the criterion for magnetization switching derived by Slonczewski and Sun [{| in their very special 
case is equivalent to the instability of a steady state in our approach. 

In another special case when we have only out-of-plane torque T±, the switching criteria are clearly the 
same as in the Stoner-Wohlfarth field-switching theory since the torque T± is equivalent to one arising 
from an effective field of magnitude H U Qvr (see Eq.(55)). The stability criteria are then equivalent to the 
conditions for a minimum of an energy function whose gradient gives all the effective fields. This was 
previously recognized by Heide et al. |16| . Obviously this energy does not involve 7 and the absence of 7 in 
both criteria is clearly seen when in the second criterion (60) we take the limit v — > 0, r — > 00 with the T± 
parameter vr finite. As soon as the in-plane torque parameter v is nonzero, no energy function exists and 
the stability criterion (60) involves the damping parameter 7, showing its essentially dynamic nature. It is 
interesting that even when both spin-transfer torques exist, the in-plane torque drops out of Eq.(60) in the 
strong damping limit 7 — > 00 and we return to a Stoner-Wohlfarth situation. 

In general, the system is neither in the Slonczewski-Sun nor Stoner-Wohlfarth limit and our general stability 
analysis based on the criteria (59) and (60) is required. We shall first apply it to discuss the stability of 
the steady-state paths shown in Fig. 3 which correspond to the microscopic parameters (b) of Table I. In 
this case, the reduced parameters of the present section, which reproduce accurately the microscopically 
determined curves in Fig. 3, are h p = and r = 1. The torque ratio r = 1 is clear from the curves for case 
(b) in Fig. 4. The corresponding steady-state paths are shown in Fig. 10 (a). In Fig. 10 (b) we plot the 
steady-state paths for the parameter set (d) of Table I, in which case r = — 1 (see Fig. 4). In Fig. 3 we 
plotted —Vf, on the bias axis, where VJ, is in volts, in order to compare with the present reduced bias v which 
is proportional to eVt and e is negative. (Vf, and v have the same sign when Tji < and opposite sign when 
T11 > 0.) We have chosen for presentation purposes a value 7 = 0.05 of the damping parameter which is 




FIG. 10: Reduced bias v required to stabilize the switching magnet moment at an angle a on the universal path for r = 1 (a) 
and r = — 1 (b) and h v = 0. Bold lines correspond to stable steady states. 

somewhat larger than that suggested in We recall that Figs. 3 and 10 correspond to the situation when 
the moment of the polarizing magnet is at 9 = 2 radians. We first discuss Fig. 10 (a). Initially with zero 
bias v the moment of the switching magnet is in the direction of the uniaxial anisotropy field, i.e., ao = 0. 
As v increases positively, ao increases until the stability criterion (60) ceases to be satisfied. In Fig. 10 stable 
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steady states are indicated by heavy lines and unstable ones by thin lines. At this value of v, the point 
A in Fig. 10 (a), the system seeks another steady state which is stable at the same bias and hence jumps 
to B. It is assumed that for finite 7 the system will home in on the stable state. Sun showed how this 
happens dynamically for the special case of 9 = 0. On further increase of bias the system proceeds to C 
and, on reducing the bias to zero, it moves to D where a — n. The current-induced switching process is 
thus completed with the magnetization switched between the two stable zero-bias orientations (a = and 
7r) in the uniaxial anisotropy field. To reverse the process the bias is reversed and the system proceeds to E 
where it becomes unstable and jumps to F. Finally a further negative increase of v takes the system to G 
and, on reducing the bias to zero, we return to O. If the resistance is calculated for each steady state, using 
our Keldysh formalism (or equivalently the Kubo formula) for charge current, the corresponding hysteresis 
loop of resistance versus bias can be plotted. We shall do this for later examples. 

In the above example, the in-plane and out-of-plane torques are of equal strength and of the same sign 
(r = 1). The instabilities at A and E are governed by the dynamical criterion (60). In the Stoner-Wohlfarth- 
like case discussed earlier the system remains stable, as bias is increased, up to the maximum at A'. The 
first criterion (59) determines the instability at this point and similarly at E'. 

Fig. 10b shows the situation for the less usual case of negative r, in particular r = — 1 corresponding to 
parameter set (d). Starting at ctQ = 0, with bias increasing from zero, an instability occurs at A as before. 
However, there are now two possible stable steady states to which the system might jump, labelled by B' 
and B. We cannot say to which point the jump occurs without following the detailed dynamics of the system 
with time-dependent solutions of the Landau-Lifshitz equation. If the system jumps to B' further increase 
of bias leads towards G where the moment of the switching magnet approaches alignment with that of the 
polarizing magnet. However the bias is varied now, through positive and negative values, the system remains 
on the stable steady state branch GG' and no switching to the point D can occur. If, however, the system 
jumps from A to B, the bias can be reduced to zero at D and a hysteresis loop can be completed via E and 
F. On the other hand if on jumping to B the bias is further increased, to reach the state C, a jump will occur 
to C and the system is again trapped on the branch GG'. 

The next example we consider employs a mapping of the fully realistic microscopic torques for Co/Cu(lll) 
shown in Fig. 7b onto the macroscopic model. In this case the switching magnet consists of two atomic 
planes of cobalt (Af = 2). We recall that the nonmagnetic spacer consists of 20 atomic planes of Cu. A 
new feature of this example is that we now introduce a strong easy-plane anisotropy with h p = 100. If we 
take H u0 = 1.86 x 10 9 sec~ 1 corresponding to an uniaxial anisotropy field of about 0.01T, this value of h p 
corresponds to the shape anisotropy for a magnetization of 1.6 x l0 6 A/m, similar to that of Co 0. We take 
again 9 = 2 radians and a realistic value of 7 = 0.01, which is in line with the value quoted by Sun The 
value of r which gives a reasonable fit to the microscopic torques shown in Fig. 7b is r = 0.65. We find that 
the strong easy-plane anisotropy forces the switching magnet moment to rotate in the (x, z) plane, which 
means that the universal paths in the (a, <fi) plane are almost straight lines with <f> = or ir. The plot of 
reduced bias v against a is shown in Fig. 11 and again the heavy portions of the curves indicate stable steady 
states. The multiple loops of steady states in Fig. 11a are a new feature appearing for h p 0. However, in 
this case, they are all unstable. The important parts of the curves are shown on a larger scale in Fig. lib. 
This clearly resembles Fig. 10a and instabilities now occur at the extremal points A' and E' instead of A 
and E as in Fig. 10a. In general, the point A lies further up the curve towards A' the larger the product jh p . 
This follows from Eq.(60) as long as the easy-plane anisotropy is strong enough for | sin^g sinao |< l/\/3 to 
be satisfied. This stabilizing effect of easy-plane anisotropy has the unwelcome consequence that the critical 
bias (current) for switching is strongly increased by such anisotropy. The corresponding hysteresis loop of 
resistance versus bias for 9 = 2 radians is shown in Fig. 12a. We have also transferred the key points from 
Fig. lib. In Fig. 12b we show the hysteresis loop for 9 = 3 radians, which is close to 9 = ir assumed in 
previous treatments Q, ||. It is rather interesting that the critical bias for switching is « 0.2mV both for 

9 = 2 and 9 = 3 radians. When this bias is converted to the current density using the calculated ballistic 
resistance of the junction, we find that the critical current for switching is ~ 10 7 A/cm 2 , which is in very 
good agreement with experiment Q. However, there is a qualitative difference between the cases 9 = 2 
and 9 = 3 radians. For 9 = 2 radians, switching is determined by the instability condition (59) which is 
independent of the damping parameter 7. This means that switching is of the Stoner-Wohlfarth type. On 
the other hand, for 9 = 3 radians, we find that the instability is determined by Eq.(60), which means that 
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FIG. 11: Reduced bias v required to stabilize the Co switching magnet moment at an angle a, with 9 = 2, and realistic 
anisotropy and damping parameters for Co/Cu/Co(lll) with M = 2. The value <f> « is obtained for v > and <j> it is 
obtained for u < 0. 
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FIG. 12: Resistance of the Co/Cu/Co(lll) junction as a function of applied bias with M = 2 monolayers of Co in the switching 
magnet, (a) is for 9 = 2 radians and (b) is for 8 = 3 radians. 

switching is of the Slonczcwski-Sun type. 

The last example we consider employs again a mapping of the fully realistic microscopic torques for 
Co/Cu(lll) shown in Fig. 7a onto the macroscopic model. In this case the number of atomic planes in the 
switching magnet is M = 1. We use the same values of h p , 7, H u o, and 9 as in the previous example. The 
best fit to the microscopic torques in Fig. 7a gives r = —1.0. This example, besides being again a realistic 
one, introduces the feature of negative r which we met in the single-orbital model with parameter set (d) 
(see Figs. 4 and 10b). The plot of reduced bias v against a is shown in Fig. 13 with stable steady states 
indicated as before. The low-bias part of the curves is plotted in Fig. 13b on a larger scale. The fact that 
| r | is comparable for Figs. 11 and 13 (of the order of 1) but the sign is changed, leads to all the bias curves 
being essentially only reflected in the a axis. However, the change of sign of r is seen to have a dramatic 
effect on the stability of the steady states. The most striking effect is the creation of a 'bias-gap' with no 
stable steady states for values of reduced bias | v | between about 0.7 and 10. For r = — 1 the bias-gap only 
exists in the presence of easy-plane anisotropy. In this connection we may compare Fig. 13, with parameters 
r = -1, hp = 100, 7 = 0.01, 9 = 2 with Fig.lOb corrcponding to r = -1, h p = 0, 7 = 0.05, 9 = 2. (The 
larger value of 7 for Fig.lOb is not important; it was used to push the point of instability A to larger bias 
and thus clarify the figure). Clearly there is no bias-gap in Fig. 10b with h p = 0. Another effect of large h p 
is to push the point of instability G in Fig. 13b to much larger bias than the corresponding point in Fig.lOb, 
even with a smaller damping parameter 7. In fact, for the particular values of 7 and h p used for Fig. 13, G 
lies at a larger bias than E'. The resultant hysteresis loop, shown in Fig. 14a, is thus executed in the same 
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FIG. 13: Reduced bias v required to stabilize the Co switching magnet moment at an angle a, with 9 = 2, and realistic 
anisotropy and damping parameters for Co/Cu/Co(lll) with M = 1. The value <j> f» is obtained for v < and (j> f» 7T is 
obtained for v > 0. 

sense as that shown in Fig. 12a. However in Fig. 14b, corresponding to 6 = 3 instead of 2, the sense is 
reversed. In Fig. 12 as we change bias from B to C or from F to G we are achieving saturation by aligning 




FIG. 14: Resistance of the Co/Cu/Co(lll) junction as a function of applied current, with M = 1 monolayer of Co in the 
switching magnet, (a) is for 9 = 2 radians and (b) is for 8 = 3 radians. 

the switching magnet parallel or antiparallel to the polarizing magnet. However, in Fig. 14 as we increase 
bias from F the steady state becomes unstable at the point G, for a critical bias, but there is no stable state 
for the system to jump to with further increase of bias (see Fig. 13b). To emphasize this point, the points G 
and C are starred in Fig. 14a. Thus for bias larger than this critical one the system cannot home in on any 
stable state and the moment of the switching magnet remains perpetually in a time-dependent state. For 
much larger bias the system can home in onto the stable multiple loop states shown in Fig. 13a. Thus there 
is a range of bias where only the time-dependent state is possible. The bottom of the gap occurs at the bias 
point G in Fig. 14a. If the bias is then reduced below this value the system will home in on a stable steady 
state and the hysteresis loop can be completed. 

We investigated the critical negative value of r at which the bias gap appears. For the parameters used 
above, the gap is not present for r = —0.05 but is already well established for r = —0.1. Thus when in-plane 
and out-of-plane spin-transfer torques have opposite sign, and easy plane anisotropy is large, only a small out- 
of-plane torque is required to produce this unusual bias gap behavior. Since out-of-plane torque corresponds 
to an effective field, we believe that this behavior is closely related to the time-dependent motion of the 
moment of the switching magnet which is observed in a sufficiently large applied magnetic field 0. This 
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alternative mechanism for time-dependent motion of the switching magnet moment is under investigation. 

Even in the absence of easy plane anisotropy, but with large negative r, we have found a critical bias 
above which only time-dependent solutions exist. However, for small negative r (| r 1) normal solutions 
can occur (no bias gap). For intermediate values of r, as discussed for r = — 1 with reference to Fig. 10b, 
switching may or may not occur. Furthermore, in the case of pure in-plane torque [r = 0) and no easy 
plane anisotropy we find straightforward switching of Slonczewski-Sun type for | cos# |> 1/3, while for 
| cos 9 |< 1/3 a more complicated hysteresis loop is found. 

CONCLUSIONS 

Our principal result is that spin-transfer torques responsible for current-induced switching of magnetiza- 
tion can be calculated quantitatively for real systems such as Co/Cu/Co junction in the ballistic regime 
using nonequilibrium Keldysh formalism. In fact, we argue that the spin-transfer torque can be calculated 
selfconsistently from first principles only in a steady state (switching magnet magnetization does not move), 
and this is precisely what the Keldysh formalism is designed for. In the small-bias (linear-response) regime 
higher-order (many-body) effects can be neglected and our results for the spin-transfer torque are, there- 
fore, quite rigorous. Keldysh formalism provides an explicit formula for the local spin current between any 
two atomic planes of the junction in terms of one-electron surface Green functions for the cut junction. 
The surface Green functions are readily available and we calculate them using a tight-binding Hamiltonian 
with parameters determined from a fit to an ab initio band structure. With the exception of Slonczewski's 
parabolic band calculation |l|, our Keldysh formulation is the only theory that yields the local spin current 
taking into into account rigorously contributions from all the parts of the junction. As the following ar- 
gument demonstrates, previous theories |l5j . which consider only scattering of spin-polarized electrons 
incident from the spacer on the spacer/switching magnet interface are incomplete. When the particle current 
flows from the polarizing magnet toward the switching magnet it is clear that a spin-transfer torque acts 
on the switching magnet. However, when the polarity of the applied bias is reversed, the current incident 
from the right lead on the switching magnet is unpolarized and, therefore, has no effect on it. On the other 
hand, it is well known experimentally |3J that changing the polarity of the bias reverses the direction of the 
spin-transfer torque (the magnitude remains the same). This cannot be explained within a theory that treats 
only the spacer/switching magnet interface. The ingredient that is missing is strong reflection of electrons 
from the polarizing magnet which results in a spin polarization of the reflected electrons. Clearly only such 
reflected spin current flowing in the direction opposite to that of the particle current can exert torque on 
the switching magnet. It follows that multiple reflections of electrons from the polarizing and switching 
magnet are an essential feature of the problem. They are treated rigorously to all orders in our theory. The 
fundamental experimental fact that the spin-transfer torque acting on the switching magnet is proportional 
to the applied bias is obtained naturally in our theory since the spin current anywhere in the junction, given 
by Eq.(25), is proportional to the difference between the Fermi functions for the left and right halves of the 
cut junction, i.e., proportional to the bias in the low-bias (linear-response) limit. These arguments indicate 
that selfconsistent treatment of the whole junction is crucial for correct understanding of current-induced 
switching of magnetization. 

The spin-transfer torque calculated from our Keldysh formalism has two components, one with the torque 
vector T|| in the plane containing the magnetizations of the two magnetic layers ('in-plane' torque) and 
another with the torque vector T± perpendicular to this plane ('out-of-plane' torque). It is generally believed 
that the effective field-like component T± is always small. We find that this is not the case and our 
calculations show that, in general, both the in-plane and out-of-plane components tend to finite values 
independent of the spacer thickness in the limit of a thick spacer. However, it is true that T± is strictly 
zero in the limit of an infinite exchange splitting between the majority and minority-spin bands in both 
ferromagnets, and this is the case considered initially by Slonczewski jl|. In the realistic case of a finite 
exchange splitting, T± is nonzero and can be comparable with Tji . The only other general case when T± 
can be small occurs for a junction with reflection symmetry about a plane at the center of the spacer. 
Hence to observe an effect of T± one needs to break the reflection symmetry of the junction. For a junction 
with polarizing and switching magnets made of the same material, this is achieved by making the switching 
magnet thinner than the polarizing magnet, and a strongest effect is found for a switching magnet only a few 
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atomic planes thick. Our calculations show that T± and Xj| are comparable for a Co/Cu/Co(lll) junction 
when the switching Co layer is one or two atomic planes thick. Nevertheless, for a good epitaxial junction 
(ballistic limit), we find that T± is 27% of Tji even for a switching Co magnet as thick as ten atomic planes. 
An alternative way to break the symmetry is to use a junction with polarizing and switching magnets made 
of different materials. 

Another result we wish to highlight is that, depending on material parameters of the junction, the relative 
sign of Tj_ and Tji can be negative as well as positive. For example, T±/T» < for Co/Cu/Co(lll) with 
a switching Co magnet of one atomic plane and T±/Tu > for two atomic planes of Co. The negative 
sign of the ratio T±/Tn has a profound effect on the stability of steady states and, hence, on the nature of 
current-induced switching. 

Finally, to determine the critical currents for switching and to investigate the effect of T± , we have used 
the microscopically calculated spin-transfer torques as an input into the phenomenological Landau-Lifshitz 
equation with Gilbert damping. Our general philosophy is that all steady states can be calculated from 
first principles and loss of their stability, determined from the Landau-Lifshitz equation, corresponds to 
switching. This holds provided there is another stable steady state at the same current density the system 
can switch into. We showed that our criterion for instability of the steady state leads to the same critical 
current for switching as that obtained earlier by Sun [f| in the special case of T± = and for the initial angle 
9 between the polarizing and switching magnet moments equal to or tt. However, we find that qualitatively 
different switching scenarios can occur when Tj_/Tm ^ 0, 9 ^ 0, and in the presence of an easy-plane (shape) 
anisotropy. In particular, when the easy-plane anisotropy is strong, even a relatively small T± (5-10 % of 
Tii) has a strong effect on switching. In the absence of an applied magnetic field, we find that an ordinary 
hysteresis loop is the only possible switching scenario when Tj_/Tj| > 0. However, for T±/Tu < 0, a normal 
hysterestic switching occurs only at relatively low current densities. When the current exceeds a critical 
value, there are no stable steady states and the system thus remains permanently in a time dependent state. 
This is analogous to the observed precession of the switching magnet magnetization caused by a DC current 
in the presence of an applied magnetic field @]. In our case, the effective field- like term Tj_, which causes 
this behavior, is proportional to the DC current and, hence, complete loss of stability of the steady state 
occurs only when this term is large enough, i.e., when the DC current is above a critical value. 

Our calculations for Co/Cu/Co(lfl) show that the critical current for switching in the hysteretic regime 
is w I0 7 A/cm 2 both for Co switching magnets of one and two atomic planes. This is in good agreement with 
experiment ;3J ■ We recall that the critical current for switching is obtained using the spin transfer torques for 
a fully realistic Co/Cu/Co(lll) junction and assuming a uniaxial anisotropy of 0.01T and Gilbert damping 
7 = 0.0 1. This is in line with the values of the uniaxial anisotropy and Gilbert damping quoted by Sun 

We conclude by stressing that all the specific results we have obtained are strictly valid for a perfect 
epitaxial junction, i.e., in the ballistic limit. However, the Keldysh formalism we have described is valid also 
in the diffusive limit. Generalization to the diffusive limit is, in principle, straightforward. For example, 
one could introduce random impurities in the lateral supercell geometry, determine the one-electron surface 
Green functions in this geometry and then perform configuration averaging of the spin current. 
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